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(N . ABSTRACT 
t-H ■ 

Context. The flux of gamma rays with energies greater than 100 MeV is dominated by diffuse emission coming from cosmic-rays (CRs) illumi- 
qq , nating the interstellar medium (ISM) of our Galaxy through the processes of Bremsstrahlung, pion production and decay, and inverse-Compton 
■ scattering. The study of this diffuse emission provides insight into the origin and transport of cosmic rays. 

' Aims. We searched for gamma-ray emission from the Small Magellanic Cloud (SMC) in order to derive constraints on the cosmic-ray population 
I and transport in an external system with properties different from the Milky Way. 
• . , Methods. We analysed the first 17 months of continuous all-sky observations by the Large Area Telescope (LAT) of the Fermi mission to deter- 
^ ■ mine the spatial distribution, flux and spectrum of the gamma-ray emission from the SMC. We also used past radio synchrotron observations of 
'"j ' the SMC to study the population of CR electrons specifically. 

\ Results. We obtained the first detection of the SMC in high-energy gamma rays, with an integrated > 100 MeV flux of (3.7 ±0.7) x 10~ 8 phcirT 2 s -1 , 
with additional systematic uncertainty of < 16%. The emission is steady and from an extended source ~ 3° in size. It is not clearly correlated with 
the distribution of massive stars or neutral gas, nor with known pulsars or supernova remnants, but a certain correlation with supergiant shells is 
observed. 

Conclusions. The observed flux implies an upper limit on the average CR nuclei density in the SMC of ~ 15% of the value measured locally in the 
Milky Way. The population of high-energy pulsars of the SMC may account for a substantial fraction of the gamma-ray flux, which would make 
the inferred CR nuclei density even lower. The average density of CR electrons derived from radio synchrotron observations is consistent with the 
same reduction factor but the uncertainties are large. From our current knowledge of the SMC, such a low CR density does not seem to be due to 
a lower rate of CR injection and rather indicates a smaller CR confinement volume characteristic size. 

Key words. Acceleration of particles - Cosmic rays - Magellanic Clouds - Gamma rays: observations 



1. Introduction 

Gamma rays with energies greater than 100 MeV are the ob- 
servable manifestations of various populations of non-thermal 



particles. The sky in this energy range was unveiled by SAS-2 
and COS-B in the 1970s, then mapped by CGRO/EGRET in 
the 1990s, and is now surveyed by AGILE since 2007 and 
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fermi/LAT since 2008. At low Galactic latitudes, many super- 
nova remnants (SNRs) and pulsars have revealed themselves as 
accelerators of particles up to TeV energies, but a lot of localised 
sources are still formally unidentified. Most of the emission, 
however, is of diffuse nature and arises from the interaction of 
Galactic cosmic-rays (CRs) with interstellar matter and radia- 
tion. The pervasiveness of this emission, which accounts for 80- 
90% of the flux above lOOMeV, indicates that CRs diffuse far 
from their sources filling the entire Galaxy. 

Since their discovery at the beginning of the 20th century, it 
has been realised that CRs are not simply a side-product of the 
most energetic phenomena of the Universe, but actually form an 
essential component of the Milky Way (MW). Galactic CRs af- 
fect the physical state of the interstellar medi um (ISM) by heat- 
ing/ionising its atomic and molecular phases dFerrier e2001) and 
are thought to alter the spectrum of interstellar turbulence as 
they get reaccel erated by draining e nergy from magnetohydro- 
dynamic waves dPtuskin et alj|2006l) . Ionisation by low-energy 
CRs is important for the chemical processes operating inside 
dense molecular clouds, where Lyman photons cannot pene- 
trate (llndriolo et al.ll200"9l) . Because of their high energy density 
(~ 1 eV cm 3 ) and corresponding pressure, C Rs also contr ibute 
to supporting interstellar gas against gravity dFerrier il998h . In 
addition, CRs are essential for the Parker dynamo mechanism, 
which may be the source of the large-scale ma gnetic fields in 
galaxies like the Milky Way dHanasz et al.| [2004. 200§). 

Direct measurements of the CR composition and spectrum in 
the interplanetary medium, together with imaging of the Galactic 
diffuse emission in the radio and MeV-GeV domain, have pro- 
vided a wealth of information on Galactic CRs. Combined with 
models of increasing sophisticati on, from simple leaky box to 
analytic diffusion models (e.g. lAharonian & Atovanl l2000b to 
more elaborate nu merical calculations (such as GALPROP, see 
Strong et al. 2004), these data permitted deriving important con- 
straints on the transport of CRs inside our Galaxy: the time spent 
in the Galactic volume, the amount of matter traversed, the ex- 
tent of spatial diffusion, a nd the role of reac celeration and con- 
vection (see the review bv lStrong et al. 1120071) . 

Resolved gamma-ray observations of externally-viewed 
galaxies provide the possibility to make further progress un- 
derstanding the origin and propagation of CRs. The Magellanic 
Clouds offer a great opportunity in this regard due their prox- 
imity and their significantly contrasting geometries and physi- 
cal conditions (gas, stars, magnetic field, ...) compared to the 
Milky Way. The Large Magellanic Cloud (LMC) was observed 
by EGRET with a measured flux >100MeV of (1.9 +0.4) x 
10~ 7 phcrrT 2 s _1 that was found to be consistent with a model 
in which CRs interact with the ISM of the galaxy, under the as- 
sumption that the CRs, magnetic pressure, and ther mal gas are in 
dynamical balance with the g ravitational attraction (Fichte l et alj 
1199 It ISreekumar etaDll992h . While it was clear that the LMC 
was an extended object, the sensitivity of EGRET precluded 
resolving any features. The LMC has now been resolved by 
Fermi/LAT. The higher sensitivity and angular resolution of 
the instrument revealed a good correlation between gamma-ray 
emission and massive star forming regions, thereby providing 
evidence for CR acceleration in these regions and suggesting 
a relat ively small diffusion length for GeV CRs (lAbdo et al.l 
l2009hl) . In contrast, the Small Magellanic Cloud (SMC) was 
not detected by EGRET but the upper-limit thus derived al- 
lowed to dismiss a metaga lactic origin for the bu lk of CRs 
with energies 10 15 -10 16 eV (ISreekumar et al.lll993h . while the 
most energetic particles are thought to come from extragalactic 
sources lying within the Greisen-Zatsepin-Kuz'min (GZK) hori- 



zon dThe Pierre Auger Collaboration et alJl2007ll2008il2009h . 

In this paper, we report the detection and resolving of the 
SMC in high-energy gamma-rays using Fermi/LAT. In the fol- 
lowing, we describe the various analyses performed on the 
Fermi/LAT data in order to derive the morphology, the flux and 
the spectrum of the SMC in gamma-rays. We then discuss plau- 
sible origins for the detected emission before considering the 
consequences of our result on the CR population and transport 
in the SMC. 



2. Observations 

2.1. Data preparation 

The characteristics and pe rformance of t h e LAT aboard Fermi 
are described in detail by lAtwood et al.l C2009). We used the 
LAT Science Tools package, which is available from the Fermi 
Science Support Center, with the P6_V3 post-launch instrument 
response functions (IRFs). These take into account pile-up and 
accidental coincidence effects in the detector subsystems that 
were not considered in the definition of the pre-launch IRFs. 

For the analysis, we selected all events within a square 
region-of-interest (ROI) of size 20° x 20° centred on {I, b) = 
(302.8°, -44.3°) and aligned in Galactic coordinated The data 
used in this work amount to 504 days of sky survey observa- 
tions over the period August 8th 2008 - December 25th 2009 
during which a total exposure of ~ 4.5 x 10 10 cm 2 s (at 1 GeV) 
was obtained for the SMC. Events of the Diffuse class (which 
are the least contam inated by the cosmic -ray background, see 
lAtwood et all 2009) and coming from zenith angles < 105° (to 
greatly reduce the c ontribution by Earth albedo gamma rays, see 
lAbdo et all 12009 gl) were used. To further reduce the effect of 
Earth albedo, the time intervals when the Earth was apprecia- 
bly within the field of view (specifically, when the centre of the 
field of view was more than 52° from the zenith) were excluded 
from this analysis. We further restricted the analysis to photon 
energies above 200 MeV because the effective area for Diffuse 
class events changes rapidly at lower energies, and up to 20 GeV 
because the number of counts is quite low beyond this limit. 

2.2. Background modelling 

At the mid-latitude position of the SMC (b ~ -45°), the gamma- 
ray background is a combination of extragalactic and Galactic 
diffuse emissions. The extragalactic component comprises re- 
solved sources, which often can be associated with known 
blazars and other extragalactic source classes, and a diffuse com- 
ponent, which is attributed to unresolved sources and eve ntually 
intrinsically diffuse pr ocesses (lAbdo et al. 2009ft l2010bh . 

The 1FGL catalog dAbdo et al.l2010al) has been derived from 

1 1 months of survey data (comparable to the data volume used 
in this analysis) and contains 622 sources for latitudes \b\ > 30°, 
corresponding to a source density of 99 sources per steradian. 
Our ROI covers a solid angle of 0.12 sr, thus we expect about 

12 resolved background sources in our field. Among those, 1-2 
should spatially overlap with the SMC for an assumed angular 
size ~ 6 - 8° for the galaxy. 

Within the ROI but outside the boundaries of the SMC, we 
found a total of 7 significant point sources, in rough agreement 
with expectations. One of these is the globular cluster 47Tuc, 



1 All the maps presented in this article are in Galactic coordinates. To 
ease the comparison with other studies of the SMC, the same maps are 
provided in equatorial coordinates as online material. 
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Fig. 1. 200MeV-20GeV counts map of the region of interest 
centred on the SMC, after smoothing with a 2D Gaussian kernel 
with cr = 0.4°. The positions of the background point sources 
are marked by white circles, with circle size indicating position 
uncertainty 



which was found to be a strong source of gamma-ray s, presum- 
ably d ue to its population of millisecond pulsars (lAbdo et al.l 
2009d). Although we found no counterparts for the other 1FGL 
sources in the CRA TES catalogue of flat-spectrum ra dio sources 
(iHealev et al.l2 007) using the procedure described in lAbdo et al.l 
(2009b), they are likely to be background blazars due to their 
high latitude positions. Searching for background sources within 
the SMC boundaries is more difficult due to possible confusion 
with point sources in the galaxy itself or with local emission 
maxima in its diffuse emission. However, background blazars 
may reveal themselves by their flaring activity. We therefore 
searched for excess emission above the time-averaged level on 
a monthly basis but found no indication for a flare coming from 
the direction of the SMC (see Sect.HTTT). 

We modelled the celestial background emission within the 
ROI by fitting components for the diffuse Galactic and extra- 
galactic backgrounds, 47Tuc, and the 6 background objects from 
the 1FGL catalogue. 

The LAT standard background model gll_iem_v02 for the 
Galactic diffuse emission was found to be inadequate for the 
specific SMC region. Subtracting this fitted model (and the 
other background components described below) from the data 
left strong positive residuals at low latitudes (b ~ -35°) that 
were found to be correlated with other tracers of the interstellar 
gas, su ch as the E(B-V) reddening map built by [Schlegel et al.l 
(119981) from COBE/DIRBE and IRAS/ISSA submillimeter ob- 
servations of dust emission (the suita bility of E(B-V ) as gas 
tracer was illustrated for instance in lAbdo et al.ll2009ib . A dedi- 
cated model for the Galactic diffuse emission in the SMC region, 
based on this E(B-V) map and an inverse-Compton (IC) inten- 
sity map computed using the GAL PROP code for CR propaga- 
tion in the Galaxy ( Stron g et al.l2000h and its model for the inter- 
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Fig. 2. 200MeV-20GeV residual counts maps after subtraction 
of the fitted celestial background model and a smoothing with a 
2D Gaussian kernel with cr = 0.4°. Upper plot corresponds to the 
full region of interest and lower plot is a zoom on the SMC. The 
positions of the background point sources are marked by white 
circles in the top panel, with circle size indicating position un- 
certainty, and the SMC is traced by H i column density contours 
in the bottom panel. 



stellar radiation field dPorter et al.ll2008l) . turned out to improve 
the situation and yielded a better likelihood. The E(B-V) map 
is converted to hydr ogen column densi ty map using the best-fit 
relation obtained by iGrenier et al.l d2005l) . and a 5° x 3° region 
around the SM C is excluded an d filled by the inpainting method 
described bv lElad et alJ (120051) . As such, it is intended to trace 
the emission due to CRs interacting with interstellar matter in 
our Galaxy, and it is fitted assuming a power-law spectral shape 
with normalisation and spectral index as free parameters. The IC 
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mapQ gives the emission arising from CRs interacting with inter- 
stellar radiation and is fitted with the spectral shape as calculated 
with GALPROP. 

An isotropic background component was fitted assuming a 
power-law spectral shape with normalisation and spectral index 
as free parameters. This component is designed to account for 
the diffuse/unresolved extragalactic emission, along with resid- 
ual instrumental background. 

The 6 background 1FGL objects are modelled as point 
sources at the positions given in the catalog and having power- 
law spectral shapes with normalisation and spectral indices as 
free parameters. Last, the globular cluster 47Tuc lying close to 
the SMC (in the plane of the sky) was modelled as a point source 
at its SIMBAD position with an exponentially cutoff power-law 
spectral shape with normalisation, spectral index, and cutoff en- 
er gy as free p arameters (in agreement with the dedicated study 
of lAbdo et al. 2009d). The characteristics of our complete celes- 
tial background model are listed in Table [1] 

3. Emission morphology 

To investigate the spatial distribution of gamma-ray emission to- 
ward the SMC, we first binned all photons into a counts map of 
size 20° x 20° centred on (/, b) = (302.8°, -44.3°) and aligned 
in Galactic coordinates. Figure Q] shows the counts map for the 
full ROI, with a pixel size of 0.05° x 0.05° and a smoothing of 
the statistical fluctuations with a 2D Gaussian kernel of width 
<x = 0.4°. The region is clearly dominated by 47Tuc close to the 
centre, and by the Galactic diffuse emission at low latitudes. 

The top panel of Fig.|2]shows the counts map for the full ROI 
after subtraction of the background model (by a binned maxi- 
mum likelihood analysis with 10 logarithmically-spaced energy 
bins), with a pixel size of 0.05° x 0.05° and a smoothing of 
the statistical fluctuations with a 2D Gaussian kernel of width 
cr = 0.4°. The map suggests that the background is properly re- 
moved by our treatment. In particular, the intensity gradient due 
to the diffuse Galactic emission increasing as we move toward 
the Galactic plane as well as the 6 background 1FGL objects and 
47Tuc have disappeared. The only remaining feature is a clear 
excess at the position of the SMC. We emphasise here that, what- 
ever the Galactic diffuse model used (the standard gll_iem_v®2 
model or our dedicated one), residual emission at the position of 
the SMC was observed and the excess counts were distributed in 
a very similar way. In addition, the characteristics of this resid- 
ual emission were found to be insensitive to the inpainting or 
interpolation method used in the preparation of the E(B-V) map 
for the dedicated model. 

The bottom panel of Fig.[2]is a zoom on the centre of the ROI 
showing residual emission over a few degrees and spatially co- 
incide nt with the SMC, as trace d by H i column density contours 
(from Stanimirovic et alJI 1999t) . In Sect. 13.11 we show that the 
extension of this excess emission is real, and not just an effect of 
the IRE The total number of excess 200 MeV - 20 GeV photons 
above the background in a 4° x 4° square covering the residual 
excess amounts to ~ 450 counts, whereas the background in the 
same area amounts to ~ 1820 counts. Such modest statistics ob- 
viously limit the precision of the conclusions about the spatial 
distribution of the extended emission from the SMC. 

The residual counts map shows several spots with intensi- 
ties of ~40-50 counts/deg 2 , similar to the peak residual inten- 
sity associated with the SMC (~65 counts/deg 2 ). None of these 

2 From GALPROP run 54_87Xexph7S, which is available from 
"http://galprop.stanford.edu". 



spots actually is significant except for the excess at (/, b) ~ 
(296°, -36°) which is at the ~ 4<x level. In contrast, the signif- 
icance of the SMC emission is ~ llcr (see Sect. 13.11 ). mainly 
because it is significantly extended while the other hot spots of 
the residual counts map are mostly point-like. We also checked 
that adding these hot spots as point sources in our model fitting 
did not alter the characteristics derived for the SMC emission. 

3. 1 . Geometrical models 

We assessed the spatial distribution of the SMC emission using 
simple parameterized geometrical models for the gamma-ray in- 
tensity distribution. We considered point-like and 2D Gaussian 
models with free geometric parameters, and for each model we 
assumed a power-law spectrum with normalisation and spec- 
tral index also treated as free parameters. In our procedure, the 
spatial and spectral parameters of the models are adjusted us- 
ing a binned maximum likelihood analysis with spatial pixels 
of 0.1° x 0.1° and 10 logarithmically spaced energy bins cov- 
ering the energy range 200 MeV - 20 GeV. For each model of 
the SMC, we computed the so-called Test Statistic (TS) which 
is defined as twice the difference between the log-likelihood -Ci 
that is obtained by fitting the model on top of the background 
model to the data, and the log-likelihood -Co that is obtained by 
fitting the background model only, i.e. TS = 2(-Ci - -Co). Under 
the hypothesis that the background model satisfactorily explains 
our data, and if the number of counts is high enough, the TS 
follows a x], distribution with p degrees of freedom, where p is 
the number of additional free parameters in the maximization of 
-Ci with respect to those used in the maximisation of -Co (in the 
present case, p is the number of free parameters of our model of 
the SMC). If, however, the TS takes a value that is statistically 
unlikely for a x\ distribution, this means that the background 
model is not an adequate enough representation of the observa- 
tions and can b e improved by the additional SMC component 
fsee lCash|[l979l) . 

We first aimed at explaining the gamma-ray emission from 
the SMC by a combination of point sources. For this purpose 
we successively added point sources to our model and opti- 
mised at each step the location, flux, and spectral index of the 
new component by maximisation of the likelihood criterion. We 
stopped the procedure when adding a new point source did not 
improve the TS value by more than 10 (which formally corre- 
sponds to a 2<x detection for an additional model component 
with 4 degrees of freedom). From this method, we obtained a 
model of the residual emission from the SMC region consist- 
ing of 3 point sources and yielding a TS of 130.4 for a total 
number of 12 free parameters. This model will be referred to 
as 3PS hereafter. The strongest point source (PS1) is located 
at (l,b) = (301.7° + 0.2°, -44.4° + 0.2°) and alone gives a TS 
of 98.2 for 4 degrees of freedom. We note here that the posi- 
tion of PS1 approximately corresponds to the position of the 
source 1FGL J0101. 3-7257 that is associated with the SMC 
in the catalog. The second point source (PS2) is located at 
(l,b) = (303.4° + 0.3°, -43.6° + 0.3°) and increases the TS to 
119.2 for an additional 4 degrees of freedom. The third point 
source (PS3) is located at (/, b) = (302.1° + 0.6°, -45.4° + 0.5°) 
and further increases the TS to 130.4 for 4 more degrees of free- 
dom. Quoted uncertainties correspond to 2cr intervals. These 
point sources correspond to the strongest maxima in the resid- 
ual emission from the SMC and their positions may mark spe- 
cific sites or objects in the galaxy: strong concentrations of CRs 
and/or gas if gamma-rays come primarily from CR-ISM inter- 
actions, or high-energy objects like pulsars. The maxima also 
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Component 


Spatial model 


Spectral model 


Parameters 


Galactic diffuse 


E(B-V) maD from Schleeel et al. (1998) 
IC map from GALPROP 


power law 
from GALPROP 


2 
1 


Extragalactic diffuse 


Isotropic 


power law 


2 


47Tuc (1FGL J0023.9-7204) 


point source at (l,b) = (305.89°, -44.89°) 


power law + cutoff 


3 


1FGL J0028.9-7028 


point source at (l,b) = (305.66°, -46.53°) 


power law 


2 


1FGL J0101. 0-6423 


point source at (l,b) = (301.22°, -52.70°) 


power law 


2 


1FGL J0217.9-6630 


point source at (/,£>) = (290.14°, -48.42°) 


power law 


2 


1FGL J2152.4-7532 


point source at (/,£>) = (315.62°, -36.82°) 


power law 


2 


1FGL J2227.4-7804 


point source at (/, b) = (311.62°, -36.51°) 


power law 


2 


1FGL J2355.9-6613 


point source at (l,b) = (311.57°, -49.96°) 


power law 


2 



could be statistical fluctuations. 

Next, instead of using a set of point sources we consid- 
ered a 2D Gaussian shaped intensity profile to build a model 
that is more appropriate for extended and diffuse emission struc- 
tures. Due to limited statistics, we restricted ourselves to a single 
Gaussian (we will see later that the photon statistics are cur- 
rently insufficient to allow a more sophisticated model). The 
best fit 2D Gaussian has a width <x = 0.9° + 0.5°, correspond- 
ing to a characteristic angular size of ~3°, and is centred on 
(I, b) = (302.1° + 0.4°, -44.4° + 0.4°). Quoted uncertainties cor- 
respond to 2cr intervals. The TS obtained for this model, referred 
to as 2DG in the following, is 136.6, which is above the TS of the 
3PS model for fewer degrees of freedom (5 instead of 12). The 
significance of the extension, when comparing the 2DG model 
with respect to the single point-source PS 1 model, is above 6<x. 
These results suggest that the emission from the SMC is diffuse 
in nature. 

Finally, we tested the possibility that the gamma-ray emis- 
sion from the SMC is not completely diffuse and has a point- 
source component, as might be expected for a pulsar adding its 
radiation to the diffuse emission arising for instance from CR- 
ISM interactions. We took this hypothetical point source to be 
the PS 1 component of the 3PS model. Essentially, we evaluated 
how the inferred diffuse emission is impacted when the max- 
imum of the residual emission from the SMC is attributed to 
a point-like object. In this case, the best-fit 2D Gaussian was 
found to be quite similar in terms of position and width, but the 
associated flux dropped by ~20% (see Table The TS value 
for this model is 138.8, a modest increase of 2.2 compared to 
the 2DG model. This attempt to improve our representation of 
the emission from the SMC actually reveals the limits on the 
spatial information that can be derived from the current set of 
observations. Yet, it also confirms that most of the emission is 
diffuse in nature: a higher TS is achieved for fewer degrees of 
freedom when complementing the PS 1 component with an ex- 
tended component rather than with more point sources, and in 
that case -80% of the flux comes from the extended component. 

3.2. Tracer maps 

The morphology study presented in !3. 11 based on analytical in- 
tensity distributions, indicates that the emission from the SMC is 
of diffuse nature. This supports the idea that the observed emis- 



sion can be attributed to CRs interacting with the ISM, although 
at this point we cannot dismiss the scenario of an ensemble of 
weak unresolved point sources in the galaxy like pulsars. To test 
the former hypothesis further, we compared our data to spatial 
templates that trace various components of the SMC relevant to 
CR-ISM interactions. In the 200 MeV - 20GeV energy range, 
the gamma-ray emission from CR-ISM interactions is domi- 
nated by the decay of neutral pions jt° created by inelastic col- 
lisions bet ween CR nuclei and interstellar gas particles (see for 
example fStrong et aLll2004l) . As tracers of CR-ISM interactions, 
we therefore used a map of neutral atomic hydrogen, the ma- 
jor target for CR nuclei in the SMC (molecular hydrogen com- 
prises <10% of the gas mass in the SMC, see 15. U . and a map of 
the ionised hydrogen that indirectly traces massive stars, that are 
putative accelerators of CRs (through their winds and/or explo- 
sions). 

For the neutral atomic hydr ogen distribution, we used the 
map of Stanim irovic et al.l l[l999) that is based on a combination 
of data from the ATCA and Parkes observatories. The map was 
converted into hydrogen column densities under the assumption 
of optically thin emission. The conversion of H i spectral line ob- 
servations to H i column densities can actually be achieved under 
various hypotheses, optically thin or thick e mission, with var- 
ious s pin temperatures in the latter case, but Stanimirovic et al. 
( 1999) showed that the correction for self-absorption in the SMC 
is a relatively small effect. In addition, the photon statistics in 
our gamma-ray data are very likely not high enough for the 
above subtleties to have an important effect on the results. The 
adopted scheme for the preparation of the H i map proved to af- 
fect only marginally the detection significance in the analysis of 
the Fermi/LAT observations of the diffuse gamma-ray emission 
of the LMC (which has a 7 times greater flux than the SMC, see 
lAbdo et al]l2009hl) . 

For the ionised hydrog en tracer, we extrac ted from the full- 
sky composite Ha map of Finkb einerl d2003l) a circular region 
of 3° in radius around the position (l,b) = (302.8°, -44.3°). A 
foreground Galactic contribution to the observed emission is es- 
timated from the off-source Ha intensity in the nearby area out- 
side our selection, and subtracted from the circular region. We 
note that a more accurate tracer of massive stars in the SMC 
could be obtained by taking into account foreground extinction 
by dust in the Milky Way and in ternal extinction in the SMC 
(see for instance [Mao et alll2008h . but we again emphasise that 
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Table 2. Comparison of the different spatial models used to represent the emission from the SMC. 



Name 


Composition 


Position 


ATS 


Parameters 


Flux 


TS 


Significance 


3PS 


point source 1 (PS1) 
point source 2 (PS2) 
point source 3 (PS3) 


(301.7°, -44.4°) 
(303.4°, -43.6°) 
(302.1°, -45.4°) 


98.2 
21.0 
11.2 


4 
4 
4 


0.93 ±0.26 
0.44 ±0.22 
0.31 ±0.19 


130.4 


9.6 


2DG 


2D Gaussian with <x=0.9° 


(302.1°, -44.4°) 




5 


2.11 ±0.24 


136.6 


10.9 


Hi 


H I column density map 






2 


2.13 ±0.25 


120.8 


10.7 


Ha 


Ha map 






2 


1.94 ±0.24 


114.4 


10.4 


PS1+ Hi 


PS1 

H I column density map 


(301.7°, -44.4°) 


98.2 
37.0 


4 

2 


0.46 ±0.27 
1.65 ±0.42 


135.2 


10.7 


psi+ Ho- 


PS1 
Ha map 


(301.7°, -44.4°) 


98.2 
29.8 


4 

2 


0.41 ±0.18 
1.47 ±0.27 


128.0 


10.4 


ps 1+ Gaussian 


PS1 

2D Gaussian with <x=1.0° 


(301.7°, -44.4°) 
(302.3°, -44.3°) 


98.2 
40.6 


4 

5 


0.39 ±0.25 
1.71 ±0.37 


138.8 


10.4 



Note to the table: From left to right, table columns are the name of the spatial model, the description of its components, the position in Galactic 
longitudes and latitudes of the point source and 2D Gaussian models, the TS increase due to each component, the number of degrees of freedom for 
each component, the flux in each component in units of 10~ 8 phcirT 2 s for photons in the 200MeV-20GeV rang e (with statistical uncerta inties 
only), the total T S of the model, and its significance in cr. The H I column density map and the Ha map are from Stanim irovic et alj d 19991) and 
Finkbeiner (2003), respectively (see text). 



the current photon statistics of our observations does not justify 
such efforts. 

Both maps were convolved with the IRFs and fitted individ- 
ually to the data (together with the background model) using 
a binned likelihood analysis with 10 energy bins and assuming 
a power-law spectral shape with free normalisation and spec- 
tral index for each model. This led to the following results: the 
Hi map gives a TS of 120.8 while the Ha map gives a TS of 
114.4, for 2 degrees of freedom in each case. The gas distri- 
bution traces the gamma-ray emission slightly better than the 
massive stars distribution. We note here that, even under the as- 
sumption of gamma-ray emission from CR-ISM interactions, H i 
or Ha distributions may not be the best representations of the 
emission. Indeed, an intensity distribution correlated with CR 
sources, such as massive stars, would mean CRs with GeV en- 
ergies are confined relatively close to their sources in a uniform 
hydrogen density. On the other hand, an intensity distribution 
correlated with the CR targets would imply a uniform CR den- 
sity. Yet, while the CR diffusion length may indeed be short 
dAbdo et al.ll200 9h), we know that in the Milky Way both the 
gas and the CRs are not uniformly distributed and this is likely 
to be the case in the SMC as well. 

As in l3.1l for the 2D Gaussian model, we tested the possibil- 
ity that the gamma-ray emission from the SMC has a composite 
origin and includes at least one strong point source. We again 
assumed this hypothetical point source to be the PS 1 component 
of the 3PS model, and we complemented it with the tracer maps. 
In this case, the PSl+Hi model gives a TS of 135.2 while the 
PSl+Ha model gives a TS of 128.0. The conclusions actually 
are the same as in 13.11 an extended component in addition to 
PS1 is more significant than other point sources like PS2 and 
PS3, and in that case ~80% of the flux is associated with that 
extended component (see Table 0. 



3.3. Correlations 

The results obtained for the various spatial models tested to ac- 
count for the SMC emission are summarised in Table [2] The 
model that best fits the data is the 2DG model, which gives the 
highest detection significance at a level of 10.9cr. Yet, the pho- 
ton statistics of the set of observations used here are not sufficient 
to conclusively distinguish the various extended emission mor- 
phologies tested in this study. 

In Fig. [3] we plot the extent of the 2DG model, computed as 
the sum of the cr of the Gaussian and the uncertainty in the posi- 
tion of its centroid, over the smoothed gamma-ray emission from 
the SMC. Although being statistically less likely than the 2DG 
model, the 3PS model gives indications on emission maxima that 
may correspond to particular sites or objects in the galaxy. We 
therefore marked these positions and their typical 2cr uncertainty 
domains in Fig. [3] For comparison, we also show in Fig.[3]the Ha 
emission contours of the SMC and the locations of the currently- 
known pulsars and SNRs^l 

The position of the maximum of the emission is consistent 
in both the 3PS and 2DG models. Even when a hypothetical 
point source is placed at the emission maximum, the best-fit 2D 
Gaussian accounting for the remaining flux is centred at about 
this same position (see model PSl+Gaussian in Table 0. From 
Figs.|2]and|3] this place does not seem to correspond to any par- 
ticular site or object in the SMC; it is neither a strong concentra- 
tion of gas or massive stars, nor the direction to one of the pulsars 
or SNRs of the galaxy. Ho wever, it is consistent wit h the centre 
of supergiant shell 304A (Stanimirovi c et al.lll999t) . located at 
(l,b) = (301.8°, -44.5°). Although they are not highly signifi- 
cant (detection at the 2-3<x level), we note that PS2 is found near 
the peak hydrogen density and close to the centre of supergiant 



3 We used the ATNF catalogue version 1.36 fro m the web page 
"http ://www.atnf.csiro.au/research/pulsar/psrcat/" ( Mancheste r et al.1 
2005) and SNRs data come from Rosa Williams web page 
"http://www.astro.illinois.edu/projects/atlas/". 
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Fig. 4. Monthly light curve of the SMC obtained with the 2DG 
spatial model and a power-law spectral shape with fixed index. 
Also shown for comparison is the light curve of 47Tuc. 



Fig. 3. Residual counts map after subtraction of the fitted celes- 
tial background model and smoothing with a 2D Gaussian kernel 
with <x = 0.4°. Logarithmically-spaced Ho- emission contours of 
the SMC are shown, together with the locations of the currently- 
known pulsars and SNRs in the galaxy (stars and points respec- 
tively). Black circles mark the extent of the 2DG model and the 
positions of the components of the 3PS model (including in all 
cases the 2cr uncertainty in the position). 



shell 37A dStanimirovic et al.|[l999h . while the position of PS3 
is consistent within uncertainties with NGC346, the most active 
star-forming region of the SMC. 

Regarding the tracers of CR-ISM interactions, both the H i 
and Ha maps are relatively good tracers of the extended emis- 
sion from the SMC. The photon statistics of the available SMC 
observations are not sufficient to conclusively distinguish one 
tracer from the other. In the case of the LMC, the gamma- 
ray emission is found to be strongly correlated with massive 
star forming regio ns and poorly correlated with gas density 



mng 

dAbdo et al.ll2009hl) . Additional observations of the SMC with 



Fermi are needed to clarify the differences and similarities in the 
gamma-ray emission of the two Magellanic Clouds. From these, 
important information will very likely be obtained about the ori- 
gin and transport of the CRs. 

4. Emission lightcurve and spectrum 

4.1. Variability 

As noted earlier, emission from within the SMC boundaries 
could arise from background blazars and be incorrectly inter- 
preted as local emission maxima of the diffuse emission from 
the galaxy. Background blazars may however reveal themselves 
by their flaring activity and we therefore searched for any excess 
emission above the time-averaged level. 

We first considered the time variability of the integrated 
gamma-ray emission from the SMC by using the extended 2DG 
model (on top of our background model). As before, the spec- 
trum of the 2DG model was assumed to be a power law, with 



the spectral index now fixed to the average value obtained from 
the fit to the entire data set. The same was done for 47Tuc to al- 
low comparison (and check that no cross-talk occurs between 
these two neighbouring sources), while the parameters of all 
other components of our background model were fixed to the 
average values obtained from the fit to the entire data set. The 
light curves of the 2DG and 47Tuc models were then derived on 
a monthly basis. The result is shown in Fig. [4] There is no indi- 
cation for a flare coming from the direction of the SMC or from 
47Tuc (in agreeme nt with the dedicated analysis of the latter, see 
lAbdo et al.ll2009dl) . 

Since background blazars obviously would appear as point- 
like sources, we also derived the light curves for the three point 
sources of the 3PS model. In this case as well, we found no in- 
dication for flaring activity from any of the three directions. 

4.2. Spectrum and flux 

So far the analysis was done assuming that the spectrum of the 
SMC emission is well described by a power law. To determine 
the spectrum of the gamma-ray emission from the SMC without 
any assumption on the spectral shape, we fitted our background 
and 2DG models to the data independently in 6 logarithmically- 
spaced energy bins covering the energy range 200 MeV - 20 
GeV. The result is shown in Fig. [5] The spectrum is relatively 
flat (F ~ 2) over most of the energy range, with a possible max- 
imum around one GeV and a break or cutoff around a few GeV. 

To determine the integrated gamma-ray flux of the spectrum, 
we fitted both simple power-law and exponentially cutoff power- 
law spectral models of the form N(E) - k (E/Eoy r exp(-E/E c ) 
to the data by means of a binned maximum likelihood analy- 
sis over the energy range 200 MeV - 20 GeV. This analysis is 
more reliable than fitting the spectra of Fig. directly since it 
accounts for the Poissonian statistics of the data. The spectral 
parameters thus obtained are listed in Table[3] The best-fit model 
is the exponentially cutoff power law. The improvement over the 
simple power-law model is ATS =6 for one more degree of free- 
dom, which corresponds to a significance of 2.4cr only. This is 
most likely because the limited statistics at high energies prevent 
the clear detection of a cutoff in the spectrum. In the following, 
we will use the flux value associated with the exponentially cut- 
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Table 3. Best-fit spectral parameters obtained from a maximum likelihood analysis with spatial model 2DG 



Model 


Index 


Cutoff 


FlOO 


F200 


TS 


Simple power law 


9 9Q+0.I0+0.07 
-0 .10-0.06 




A tq+0.71+0.98 
t-JO_0.49-0.98 


9 | 1 +0.24+0.12 
^ -0.24-0.12 


136.6 


Exponentially cutoff power law 


l 7 £+0.22 +0.00 
" 1- -0.14-0.01 


t a+3.6+1.8 
- , -1.3-0.8 


t ££+0.74 +0.34 
-0.66-0.60 


1 nr)+0.26+0.11 
'* -0.26-0.13 


142.6 



Note to the table: From left to right, table columns are the spectral model, the spectral index, the cutoff energy (in GeV), the photon flux in 
the 100MeV-500GeV and 200MeV-20GeV range (both in 1CT 8 phcnT 2 s~'), and the TS value. Quoted errors are statistical and systematic 
respectively. 



off power-law model because it is more relevant to the physical 
sources we will consider than the simple power-law model: we 
indeed discuss CRs-ISM interactions as a likely source of the 
SMC emission and the corresponding radiation is usually domi- 
nated by a ~GeV bump due to tt° decay, and this feature is better 
represented with an exponentially cutoff power-law model (see 
Sect. 15.11 and Fig. |5j; alternatively, we consider pulsars as an- 
other plausible source of the SMC emission and the typical spec- 
trum of the latter is characterised by a cutoff at a few GeV (see 
Sect. 15.2b . We note that the relatively small difference between 
both spectral models ensures that the morphological analysis of 
the SMC emission presented in [3] and performed with an as- 
sumed simple power-law spectral shape, remains valid. 

The integrated flux of the best-fit model is (1.9 ±0.2) X 
lO^phcm^s" 1 for photons in the 200MeV-20GeV range, 
and (3.7 +0.7) x 10~ 8 phcirr 2 s_1 for photons in the 100 Me V- 
500 GeV range, all uncertainties being purely statistical. Since 
we restricted our analysis to photons > 200 Me V the extrapola- 
tion down to lOOMeV may introduce a systematic uncertainty 
if the low-energy spectrum of the source deviates from our fit- 
ted model. In addition, systematic uncertainties on the LAT in- 
strument response function translate into systematic uncertain- 
ties on the spectral parameters that are given in Table [3] The 
> 100 Me V flux is provided to facilitate comparison wi th the 
upper limit from CGRO/EGRET. ISreekumar et all dl993h gives 
an upper limit of 5.0 x 10~ 8 phcirT 2 s from the first years of 
observations of the SMC with CGRO/EGRET, which is com- 
patible with our extrapolated flux, whatever the spectral model 
used. 



5. Gamma-ray sources in the SMC 

Cosmic-ray interactions with the ISM may well be the primary 
source of gamma rays in the SMC. The resulting emission is 
expected to be intrinsically diffuse in nature and our analysis 
indeed suggests that the gamma-ray emission from the SMC 
does not originate in a small number of individual point sources. 
When trying to account for the emission from the SMC with 
a combination of point sources, we detected only three point 
sources with significance above 2cr within the boundaries of the 
galaxy. In contrast, we found a model with less degrees of free- 
dom that resulted in a larger TS value when we allowed for the 
source to be extended (model 2DG). The SMC emission is there- 
fore more likely diffuse in nature. Alternatively, it could be com- 
posed of a large number of unresolved and faint sources that can 
not be detected individually by Fermi/LAT, In the following, we 
examine these two possibilities in turn. 



5.1. CRs - ISM interactions 

In this section, we make the assumption that a large fraction of 
the gamma-ray emission from the SMC indeed originates from 
CRs interacting with the interstellar gas and radiation field. We 
first computed the average integrated > 100 Me V gamma-ray 
emissivity per hydrogen atom of the SMC using the following 
relation: 

m _ 

2 (1) 



>100MeV ( j ) >100MeV 



M, 



gas 



where <t> y is the integrated > 100 Me V photon flux from the 
SMC, M gas is the total gas content of the galaxy, d is the distance 
to the galaxy and m p is the proton mass. In a more practical form 
(where >100MeV is implied): 



q y = 8.010 



■30 



lO^phcm^s- 1 \lkpcj \10 8 M ( 



—\ 2 



M„, 



(2) 



where the emissivity q y is in phs 1 sr~' H _1 . With Eq. [TJ and 
[2] we focus on the gamma-rays produced by the interaction of 
CRs with interstellar matter. These are produced by two pro- 
cesses: decay of neutral pions n° created by inelastic collisions 
between CR nuclei and ISM gas particles, or Bremsstrahlung 
of CR electrons and positrons. In the local Galactic conditions, 
the contribution of Brems strahlung becomes small beyond a few 
100 Me V (see for instance lAbdo et alj2009eh . and we will there- 
fore consider in the following that gamma-rays from interactions 
with interstellar matter are mostly due to the hadronic compo- 
nent of CRs. The questions of the average electron spectrum and 
electron-to-proton ratio in the SMC are addressed in Sect. 16.11 
andO 

The SMC is a gas-rich system (compared to typical spiral 
galaxies, see lHunterifl997T ). The total mass of neutral atomic 
hydrogen (HI) has been determined from 21cm radio obser- 
vations with the Parkes and ATCA facilities; after correction 
for self-absorption, it w as found to amount to 4.2 xl0 8 M o 
(IStanimirovic e t al. 1999). The total mass of molecular hydrogen 
(H2) has been determined from far-infrared observations with 
the Spitzer space telescope and found to be 3.2 xlO 7 M . The ra- 
tio of molecular to atomic gas is ther efore 2-3 times lower than 
in spiral galaxies dLerov et al. 2007). We thus obtain a HI+H2 
gas mass of ~ 4.5 x 1O 8 M . This estimate, combined with the 
gamma-ray flux found in 14. 21 translates into an average gamma- 
ray emissivity > 100 MeV of (2.5 +0.5) x lO ^ph s' 1 sr' 1 H' 1 
for an adopted distance to the SMC of 61.5 kp c dCioni et al.l 
2000; iHilditch et al]l2005t iKeller & Woodll2006t) . with the un- 
certainty reflecting only the uncertainty on the gamma-ray flux. 

The average gamma-ray emissivity of the SMC is a factor of 
6-7 times smaller than the emissivity q y - (1.63 ± 0.05) x 10~ 26 
phs'sr^H 1 that has been det ermined with Ferm i/LAT for 
the local ISM of our own Galaxy dAbdo et al.ll2009eh . However, 
these quantities are not directly comparable for the reasons we 
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list below. 

First, our derivation of the average gamma-ray emissivity per 
H-atom of the SMC implicitly assumes that all gamma-rays are 
produced by n decay after hadronic interactions. Yet, we know 
from the study of the Galactic high-energy diffuse emission that 
another process of leptonic origin, the inverse-Compton scatter- 
ing, can account for a certain part of the flux in the 200 MeV 
- 20 GeV energy range. The contribution of inverse-Compton 
was subtracted prior to the determination of the local emissivity 
using predictions from the GALPROP model, but the average 
emissivity of the SMC was computed from a mix of n° decay 
and inverse-Compton gamma-rays in unknown proportions. 

Second, some of the observed emission from the SMC may 
arise from unresolved point sources like pulsars, and this would 
overestimate the flux due to CRs interacting with interstellar 
matter (see Sect. I5.21 i. In contrast, the local emissivity was de- 
termined for high-latitude regions of the Galaxy, mostly within 
1 kpc from the solar system. The observations were cleaned for 
point sources (at least for the point sources detected by the LAT 
over the first 6 months) and it seems reasonable to assume that 
the emission from the small volume under investigation was left 
relatively free of contamination by isolated objects. 

Last, the emissivity is computed from the mass of hydro- 
gen and as evident from Eq. [1] the higher the mass for a given 
gamma-ray flux, the lower the emissivity. The value we used 
for the SMC (4.5 10 8 M ) is probably a lower limit of the total 
mass of gas in the galaxy. First, we neglected the ionised gas, 
which can account fo r up to 10% (in the case of the MW, see 
ICordes & Lazioll2002l) . More important is the possibility that the 
total SMC gas mass is not completely traced by the usual meth- 
ods. This additional gas could be in the form of cold atomic gas 
being optically thick (thus not detectable at 21 cm), or molecular 
gas free of CO (thus not detectable at 2.6 mm). The latter sce- 
nario would be particularly relevant for a low-metallicity system 
such as the SMC, in which CO may be underabundant and/or 
not efficiently shielded from dissociating radiation because of a 
low dust content. In contrast, the high-latitude observations of 
the solar neighbourhood are probably not too much affected by 
this dark gas because the latter surrounds the CO clouds and 
is thus mostly confined to the Galactic pla ne. In their analysis 
of a region of the second Galactic quadrant, Abdoetal- (2009i) 
found that most o f the dark gas is con fined to Galactic latitudes 
\b\ < 20°, whereas lXbdo et alJ (f2009e) derived the local emissiv- 
ity from a region of the third Galact ic quadrant with \b\ > 22°. 
In addition, it is worth noting that lAbdo et al .] <[2009e) found 
good agreement between the emissivity inferred from the dif- 
fuse gamma-ray emission and the emissivity computed from the 
measured spectrum of CRs without resorting to any dark gas 
contribution. 

On the whole, this means that the true average emissivity for 
the SMC is very likely below the above-mentioned value. The 
average gamma-ray emissivity of the SMC should therefore be 
regarded as being at least 6-7 times smaller than the local value. 
We emphasise, however, that depending on the actual relative 
distributions of gas and CRs in the SMC, the emissivity may 
reach higher values in some places. This was clearly illustrated 
in the recent Fermi/LAT study of the LMC. The average emis- 
sivity of the LMC is found to be 2-4 times smaller than the local 
value. Yet, the gamma-ray emission from the LMC is observed 
to be poorly correlated with the gas density, and this translates 
into a gamma-ray emissivity being higher than th e local value in 
several places (see Fig. 10 in lAbdo et al] |2009h). For the SMC, 
the current photon statistics are insufficient to build an emissiv- 
ity map of the galaxy, so we must confine ourselves to the global 



average. 

In comparing the gamma-ray emissivities per H-atom of the 
local Galaxy and the SMC, one should also not forget the effect 
of the chemical composition of both the ISM and the CRs. The 
latter is usually taken into account through the so-called nuclear 
enhancement factor e^, which is a multiplicative factor applied 
to the emissivit y due to pr oton-proton reactions only. Recent cal- 
culations bv lMoril ([2009) yielded = 1.85 for CRs having the 
composition measured at Earth and interacting with an ISM of 
solar metallicity. The Z= 0.2 Z metallicity of the SMC dDufourl 
1984) corresponds to = 1.72, so we can consider that the ef- 
fect of the low-metallicity of the SMC on the gamma-ray emis- 
sivity is negligible. 

Our result, together with the above arguments, indicates that 
the average density of CR nuclei in the SMC is at least 6-7 times 
smaller than in the MW at the position of the solar system if the 
energy spectrum of both populations of CR nuclei are assumed 
similar. Yet, there is no reason for the CR spectrum not to vary 
from one galaxy to another. Indeed, the spectrum of the large- 
scale, steady-state population of CRs is very likely shaped pri- 
marily by transport, since we do not expect the diffusive shock 
acceleration mechanism to be strongly different, and transport 
of CRs may be highly dependent on characteristics of the host 
galaxy like geometry, magnetic field strength and topology, or 
level of interstellar and magnetic field turbulence. The actual 
spectrum of CR nuclei can be derived from the shape of the 
gamma-ray spectrum. However, our data do not have sufficient 
statistics to be conclusively constraining in this regard. These is- 
sues are addressed in more detail in l6.2l through modelling of the 
gamma-ray emission from CR-ISM interactions in the SMC. 

5.2. Isolated sources 

The current Fermi/LAT observations of the SMC are limited in 
their capability to rule out the possibility that a sizeable fraction 
of the diffuse gamma-ray emission is due to a large number of 
unresolved and faint sources. In this section, we estimate that 
contribution to the observed emission. 

The 1FGL catalog contains a total of 1451 sources, about 
half of which were associated with candidate sources. Most of 
the associated sources are active galactic nuclei, and the remain- 
ing as sociated sources a re mostly Galactic objects (see for in- 
stance Abdoetal . 2009b). Although there is still a large number 
of unidentified sources, many of which lie in the Galactic plane, 
the Galactic population of gamma-ray sources seems to be dom- 
inated by pulsars. In the following, we will therefore focus on 
this class of object as possible sources of gamma-ray emission 
in the SMC. 

The first Fermi/LAT catalogue of gamma-ray pulsars based 
on the first six months of LAT observations contains 46 objects, 
most o f which are normal non-recycled pulsars (Abdo et al. 
l2009al) . From distance estimates and using a flux correction fac- 
tor fn = 1, the observed fluxes were converted into gamma- 
ray luminosities > 100 MeV approximately ranging from 10 34 to 
10 36 erg s _1 for most normal non-recycled pulsars. The inferred 
luminosities suffer from large uncertainties on the distance esti- 
mates, which are manifest from the > 100% efficiencies implied 
for several objects. Another source of error is the assumption 
of a uniform beaming factor /n = 1 across the sky, which is 
not realised in many emission models and thus overestimates the 
gamma-ray luminosity in some cases. The characteristic ages for 
most normal non-recycled pulsars range from 10 4 to 10 6 yrs, and 
their typical spectrum has been found to be an exponentially cut- 
off power law with index 1.4 and cut-off energy 2.2 GeV. In the 
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following, we will assume that normal gamma-ray pulsars have 
a mean luminosity of 10 35 erg s _1 and a mean lifetime of 10 5 yrs. 
Adopting the above mentioned mean spectrum, the energy lu- 
minosity of 10 35 ergs~' translates approximately into a photon 
lu minosity of 10 38 ph s" 1 . 

ICrawford"eTaiT(l200l estimated the number of active radio 
pulsars in the SMC to 5100 ±3600, under the assumption that 
radio pulsars form at the same rate as SNII and have a mean 
lifetime of lOMyr. This would imply 51 +36 gamma-ray pul- 
sars, under the hypothesis that all young pulsars are gamma- 
ray emitters with a lifetime of 0. 1 Myr. We currently know 6 
pulsars in the SMCfl 5 discovered in rad io and 1 in X-rays 
ilManchester et all 120051: lLamb et alJl2002h . Most of these pul- 
sars have characteristic ages above 1 Myr and may therefore be 
too old for substantial gamma-ray emission. A possible excep- 
tion is JO 100-721 1, which is an Anom alous X-ray Pulsar can di- 
date with an estimated age of 6800 yrs dMcGarrv et alj|2005l) . 

A mean luminosity of 10 38 phs~' per pulsar in the SMC 
translates into a photon flux at Earth of ~ 2 10 _10 phcm -2 s _1 . 
The 51 +36 gamma-ray pulsars would therefore give rise to a 
gamma-ray flux of (1.0 + 0.7) 10~ 8 phctrT 2 s _1 , to be compared 
with the observed value of (3.7 + 0.7) 10~ 8 phctrT 2 s -1 . Taking a 
lifetime of 1 Myr instead of 0. 1 Myr for gamma-ray pulsars leads 
to the conclusion that the entire SMC emission can be accounted 
for by pulsars. Yet, these estimates suffers from large uncertain- 
ties, first in the birth rate of gamma-ray pulsars and then for the 
gamma-ray emission from a typical pulsar population. A more 
reliable calculation would require a complete population syn- 
thesis including, in particular, the time evolution of the gamma- 
ray luminosity. However, such a simulation is quite sensitive to 
the as yet unsettled magnetospheric pr ocesses governing high- 
energy emission (see lAbdo et alj|2009al and references therein). 
Nevertheless, the above calculation shows that a substantial frac- 
tion of the observed emission could be accounted for by normal 
pulsars. Adding other objects like millisecond pulsars would fur- 
ther increase the contribution of isolated sources. 

On the whole, this means that the density of CRs in the SMC 
may be even smaller than the values discussed previously, and 
hence significantly below the local value. This also implies that 
the gamma-ray emission of the SMC may be of a different com- 
position compared to the Milky Way. In the Galaxy, only 10-20% 
of the gamma-ray emission can be attributed to isolated sources 
while the rest is of diffuse nature, but in the SMC the proportions 
may be inverted. Such a scenario would require a lower mean 
SN rate per unit volume (to form gamma-ray pulsars without in- 
creasing the mean CR injection rate per unit volume) and/or a 
lower confinement of CRs (to reduce the accumulation of parti- 
cles in the ISM). We will discuss this in more detail in !6.3l 

On that issue, we note that the refinement of the emission 
morphology that will result from increased exposure may al- 
low the assessment of the relative contributions of each kind of 
source. Also, the typical spectrum of pulsars has a specific shape 
that may not be consistent with the better-defined high-energy 
spectrum that will follow from better statistics. 



4 These 6 objects actually are those listed in the ATNF catalogue re- 
ferred to previously, but there are in fact many more known pulsars in 
the SMC. Galac he et all i2008h lists about 30 X-ray pulsars that are pre- 
sumably parts of Be/X-ray binaries. Yet, these pulsars are not expected 
to be powerful gamma-ray emitters because of their long periods >10s. 



6. The CR population of the SMC 

In this section, we focus on the implications in terms of CR pop- 
ulation of the detection of the SMC in gamma-rays. We again 
assume that the entire emission is due to CR-ISM interactions 
and we derive the characteristics of the CR population of the 
SMC using a more detailed approach than the simple emissivity 
calculation performed in Sect. 15.11 In order to broaden the is- 
sue, we investigated the population of CR electrons in the SMC. 
Compared to CR nuclei, CR electrons undergo severe radiative 
energy losses. While the main loss mechanism for nuclei above 
1 GeV is escape from the host galaxy, the CR electrons lose the 
bulk of their energy through inverse-Compton scattering of the 
interstellar radiation field and synchrotron radiation on the galac- 
tic magnetic field. Hence, the most energetic CR electrons do 
not travel far from their source regions before losing most of 
their energy. By considering the conditions for the CR nuclei and 
electrons in the SMC, it may be possible to assess if the lower 
density of CRs in the SMC is due to injection or to transport. In 
addition, we discuss the aspects of injection and/or transport that 
may drive the large-scale properties of the CR population of the 
SMC and their differences with respect to the Galaxy. 

6. 1 . Constraints from synchrotron emission 

When accelerated in the large-scale magnetic field of a galaxy, 
CR electrons manifest themsel ves through diffuse synchrotron 
emission. From Longair ( 1994), the synchrotron radiation spec- 
trum of a power-law energy distribution of electrons is (with 
some changes in the units): 

J(v) = 2.344 x 10^ f(p)B* K e l ^ 217 * 10? ) V (3) 

where v is the radiation frequency in Hz, J(v) is the emissivity 
per unit volume in WirT 3 Hz -1 , B is the magnetic field in /iG, 
N e {E) = K e {E I Eq)~ p is the energy spectrum of CR electrons 
in cirT 3 GeV -1 and f(p) is a numerical factor depending on the 
slope of the particle spectrum. Then, the intensity of a given ob- 
ject is given by: 

dQ. 

F(v) = J(v) — dl (4) 
4n 

where dQ. is the solid angle subtended by the diffuse source, in 
sr, and dl is the mean depth of the emitting medium along the 
line of sight, in m. The intensity F(v) is in WirT 2 Hz 1 (which 
corresponds to 10 26 Jy). 

The radio continuum emission from galaxies usually has a 
steep spectrum at lower frequencies, due to synchrotron emis- 
sion from non-thermal particles, followed by a flatter spectrum at 
higher frequencies, due to free-free emission from therm al par- 
ticles (with a typical index of -O.ll. lHavnes et all (Il99ll) deter- 
mined a spectral index of -0.85 +0.10 for the non-thermal radio 
continuum emission of the SMC, with a flux densit>0of 34 +6 Jy 
at 1.4 GHz. 

From Eq. [3] the radio measurements give p- -2.7 +0.2 for 
the average spectral index of the CR electron distribution of the 
SMC. The normalisation of the CR electron spectrum can be ob- 
tained from the flux density at a given frequency if we know 
the magnetic field strength and the spatial distribution of the 
synchrotron-emitting medium (see Eq.[3]). 

5 The total flux density at 1 .4 GHz is 42 ±6 Jy and we subtracted 20% 
for the contribution of thermal emission (as estimated by the authors at 
1 GHz). 
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Th e large-sca l e mag netic field of the SMC was recently stud- 
ied by iMao et al.l (|2008) using optical starlight polarisation and 
rotation measures for background extragalactic radio sources. 
The authors determined that the coherent magnetic field in the 
SMC has a strength of ~ 1.7yuG and lies almost entirely in the 
plane of the sky. The SMC also seems to have a random mag- 
netic field with an estimated strength of ~ 2-3 fiG depending on 
the method used. The observed synchrotron emission from the 
SMC actually arises from the total magnetic field projected onto 
the plane of the sky B tot ai,±, and lMao et alJ (f2008) estimated this 
component to have a strength of ~ 3 yuG. 

The 3D structure of the SMC is quite controversial and this 
directly affects the interpretation of the radio synchrotron emis- 
sion. The size and appearance of the SMC in the sky is biased by 
observations of young stars and gas, which suggest a very irregu- 
lar and asymmetric shape while the old stellar population s eems 
to follow a very regular distribution (Gonid akis et alj20 09). Yet, 
for the purpose of interpreting radio synchrotron emission from 
CR electrons, young, massive stars are more relevant than old 
stars since they are putative sources of CRs. From a visual in- 
spection of the radio maps of the SMC, it appears that the emit- 
ting regions can be enclosed in a circle of 5° in diameter cor- 
responding to a solid angle of dQ. =6 10~ 3 sr. However, the ma- 
jority of the emission likely comes from half of this solid an- 
gle. The issue of the extension of the SMC along the line of 
sight is important as well. Past estimates of the depth of the 
gal axy range from a few kp c to a few lOkpc (see the review 
by [Stanimirovic et a"D l2004l section 4.2). Currently, lower val- 
ues are favoured. The deviations of pulsating red giants in the 
SMC from the mean period-luminosity relations indicate a dis - 
tribution over a distance range of 3.2 + 1.6kpc (lLah et al.i r2005). 
This resul t is reminiscent of the distance dispersion of ± 3.3 kpc 
found bv lWelch et all d 1987b from an analysis of Cepheid vari- 
ables in the SMC (and put forward as an argument that the SMC 
is not extended beyond its tidal radius and is therefore not un- 
dergoing disintegration). In the following, we will adopt a mean 
depth along the line of sight of dl =4 kpc. In our further calcu- 
lations, however, w e will also consid er larger depth such as the 
10 +6 kpc found bv lMao et alJ d2008h when deriving the line of 
sight magnetic field of the SMC from rotation measures. 

From the above data and Eqs[3]and|4] we computed the nor- 
malisation factor for the CR electron population of the SMC and 
obtained K e = 2.5 10 2 cm -3 GeV -1 . This corresponds to an 
energy density of wcu, e = 3.6 10 _3 eVcrrT 3 for particles above 
1 GeV with a spectral index of p- -2.7. For comparison, the lo- 
cal energy density of CR electrons estimated from the recent 
Fermi/LAT measurement is wcu, e - 6.0 10~ 3 eVcm " 3 for parti- 
cles ab ove 1 GeV with a spectral index of p- -3.0 ( Abd o et aD 
2009c, assuming the spectral slope measured over the 20 GeV- 
1 TeV range extends down to 1 GeV). 

The average density of CR electrons in the SMC from this 
argument is ~ 60% of the Galactic local value. Yet, when taking 
into account the uncertainties on the radio data (spectral index 
and intensity), we find that the CR density derived for the SMC 
undergo relative variations from -50% to +70%. Also, the depth 
of the galaxy may be larger than the 4 kpc we adopted and as- 
suming a depth of 10 kpc instead of 4 kpc would decrease the 
CR electron energy density by 60%. On the other hand, the size 
of the synchrotron emitting region may be smaller by 50% com- 
pared to the value used so far and this would nearly compensate 
an increased depth. 

Large uncertainties therefore affect the derivation of the CR 
electron density in the SMC from radio synchrotron observa- 
tions. At present, it is consistent with the local value and with 



a 6-7 reduction factor as observed for CR nuclei (the latter case 
implying the same electron-to-proton ratio as measured locally). 
CR electrons also manifest themselves through inverse-Compton 
gamma-ray emission, which offers an alternative to address the 
above questions. However, as discussed in the next section, the 
statistical limitations of the current high-energy observations are 
not very constraining either. 



6.2. Constraints from gamma-ray emission 

We simulated a gamma-ray spectrum for the SMC using a one- 
zone model of CR-ISM interactions that takes into account 
7T° decay following proton-proton interactions, Bremsstrahlung 
from CR electrons and inverse-Compton scattering of cosmic- 
ray electrons on optical, infrared and cosmic microwave back- 
ground photons. In this model, the various quantities involved in 
the CR-ISM interactions are represented as average values over 
a given volume or gas mass, depending of the emission process. 
For the CR populatio n, we used the proto n, electron and positron 
spectra presented in lAbdo et al. (2009e) for the local Galactic 
environment, with a global scaling factor to account for the CR 
density difference between the SMC and the local MW (see be- 
low). 

We calculated the tt° production by proton-proton interac- 
tions following the prescription of iKamae et al.l (120061) . The n° 
and Bremsstrahlung emissivities were c alculated assuming an 
average SMC metallicity of Z = 0.2 Z (lDufourlll984l) , which 
implies a nuclear enha ncement factor of - 1.72 for the n° 
emissivity (Mori 2009). These average emissivities are then mul- 
tiplied by the the total hydrogen mass of the SMC (4.5 x 10 8 M , 
see Sect. 15. U to give the total n° and Bremsstrahlung gamma-ray 
luminosity. 

The inverse Compt on component was c a lculate d from the 
method described by iBlumenthal & Gould! d 19701) using the 
CMB, optical and infrared interstellar radiation fields (ISRFs). 
The optical and infrared ISRF s were taken from the GALPROP 
model of the MW dPorter et al.ll2.008h and averaged over a cylin- 
drical volume with a 15 kpc radius and a 4kpc thickness. The 
resulting average optical and infrared ISRFs were rescaled ac- 
cording to the stellar luminosity density and observed infrared 
emis sion for the SM C (factors of -1.2 and ~1.25 respectively, 
see lBot et al.ll200 4). The inverse Compton gamma-ray luminos- 
ity was then calculated assuming that the CR electrons of SMC 
interact with the target photons in a cylinder with a 1.6 kpc ra- 
dius (radius of the SMC gamma-ray emission, see 13. U and a 
4 kpc height (depth of the SMC along the line of sight, see 16. 11 1. 

The spectral profile thus obtained is associated with the 2DG 
intensity distribution to form a complete model for the SMC 
emission. The scaling factor r c of our model with respect to the 
data, which we refer to as the cosmic -ray enhancement factor, is 
then a direct measure of the average density of CRs in the SMC 
with respect to that in the vicinity of the Earth. From a binned 
maximum likelihood analysis, we found r c = 0.14 with a statis- 
tical uncertainty of +0.01. Systematic errors due to uncertainties 
in the effective area of the instrument amount to +0.01, and an 
additional uncertainty of ±0.01 comes from the uncertainty on 
the extent of the source (the <x of the 2D Gaussian model). The 
r c factor obtained from this spectral modelling is consistent with 
the estimate made in 15. II 

From this one-zone model of CR-ISM interactions in the 
SMC, the assumption of CR spectra and CR electron-to-proton 
ratio identical to the locally observed values does not conflict 
with the gamma-ray spectrum. 
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Fig. 5. Spectrum of the SMC. The LAT data points (in black) 
come from independent fits of the 2DG plus background models 
in 6 different energy bins. The curves correspond the compo- 
nents and total of a spectral model of CR-ISM interactions in 
the SMC fitted to the data using a binned maximum likelihood 
analysis (see text). 



6.3. Injection and transport 

The average density of CR nuclei in the SMC seems to be at least 
6-7 times smaller than in our local environment, if we assume a 
similar spectrum for the particles. The energy distribution of CR 
nuclei is poorly constrained by our gamma-ray spectrum. The 
density of CR electrons in the SMC was found to be ~ 2 times 
smaller than in our local environment, but due to uncertainties 
on the radio data and other parameters, it could be down to 10 
times smaller or consistent with the local value. The specifics of 
CR injection and transport in the SMC with respect to the Milky 
Way lead to lower average CR densities and possibly also to a 
higher electron-to-proton ratio, although a similar electron-to- 
proton ratio cannot be excluded. 

Similar results were obtained in a recent study of the LMC 
based on Fermi/LAT observations. We report here the main con- 
clusions of that work for comparison. The gamma-ray emis- 
sion from the LMC is clearly of diffuse nature, with an exten- 
sion of angular size ~ 3.5°. The total flux from the LMC is 
(2.6 +0.2 ±0.4) x 10~ 7 phcm~ 2 s _1 , which is about 7 times the 
flux we obtained for the SMC. The spectrum was found to be 
well represented by a power law with an index of about -2.0 and 
a cut-off at a few GeV, quite similar to what we obtained for the 
SMC. These observations, combined with the properties of the 
LMC, imply a CR density that is 2-4 times smaller than the lo- 
cal Galactic value, possibly even less. Regarding CR electrons in 
the LMC, the authors did not explicitly estimate their density but 
they showed that the assumption of a CR electron-to-proton ra- 
tio identical to the locally observed value does not conflict with 
the gamma-ray spectrum or the radio synchrotron observations 
of the LMC. 

A lower density of CR nuclei in a steady-state system can 
result from a lower rate of CR injection and/or a reduced con- 
finement in the volume of the galaxy. To first order, the injection 
rate can be simply re lated to the rate of supe rnovae or to the 
rate of star formation dPavlidou & Fieldsll2001l) . because we do 
not expect the diffusive shock acceleration process occurring in 



supernova explosions to strongly differ from one galaxy to an- 
other. Conversely, CR confinement may be highly dependent on 
some characteristics of the galaxy like geometry, magnetic field 
strength and topology or the presence of a galactic wind and/or 
halo. In the following, we discuss those specifics of the SMC 
that may explain the lower CR densities. 

Th e star formation history derived by lHarris & Zaritskvl 
(2004) from UBVI photometry of over 6 million SMC stars show 
that the star formation rate (SFR) in the SMC over the last Gyr 
has an underlying constant value of ~0. 1 M yr _1 , with episodes 
of enhanced activity at 60 and 400 Myr where the value was 
increased by factors of 2-3. Based on far -infrared dust emis- 
sion and Ha emission, IWilke et al l (|2004) deduced a present- 
day SFR of ~0.05 M H yr 1 , in good agreement with the value of 
lHarris & Zaritskvl (|2004j) and with previous estimates. By com- 
parison, a review of the SFR estimates for the Milky Way by 
Diehl et all (120061) shows that a value of 4 ±2M yr~' can be 
adopted. The SFR in our Galaxy is therefore about 40 times 
higher than in the SMC. We note that a similar scaling is found 
for the supernova rates: the commonly adopted value for the 
Milky Way is ~3 per century, while the es timates for the SMC 
are ~0.1 per century (Tamma nn et al.l[T994l) . 

Yet, for the issue of steady-state CR density in a given sys- 
tem, it is not the absolute rate of star formation or supernovae 
that matters, but the rate by unit volume. Estimating the vol- 
umes in which star formation occurs is not an easy task, es- 
pecially for the SMC, the exact morphology of which is still 
unclear (see Sect. 16.1b . So instead of trying to compute these 
volumes, we used the interstellar radiation field density as an in- 
dicator of the star formation per unit volume. The average tem- 
perature of the dust is observed to be higher in the SMC than in 
the Milky Way, suggesting a stronger inters tellar radiation field 
(IStanimirovic et aLlEoOOt Iwiike et alj|2004l) . This could be due 
to the lower metallicity of the SMC, but may also be the con- 
sequence of a higher SFR per unit volume. If the rate of CR 
injection per unit volume actually is higher in the SMC than in 
the Galaxy, then the average lower CR densities we inferred for 
the SMC could be due to some properties of their transport in 
the ISM. 

In a basic picture, the diffusion of CRs in the ISM is gov- 
erned by two parameters: the spatial diffusion coefficient D xx 
and the confinement volume characteristic size l co „f (both pa- 
rameters being energy-dependent). The main ISM characteristic 
thought to affect D xx is the B/SB parameter (see the review by 
IStrong et alJ2007l) . From the estimate of the random and ordered 
magnetic field strengths in the SMC bv lMao etal . (2008), there 
is no reason to think that D xx in the SMC could be substantially 
different from D xx in the Galaxy. We are therefore left with the 
confinement volume size l con f as the potential explanation for 
the lower CR density of the SMC. Below, we make a first-order 
comparison of this parameter for the SMC and MW. 

Neglecting energy losses, the steady-state average CR den- 
sity in a galaxy can b e approximated by (see for instance 
iPavlidou & Fieldsll200ll) : 

N(E) = q(E)T esc (E) (5) 

where T esc is the characteristic escape time and q(E) the source 
term. For a diffusive process, the confinement volume size can 
be approximated by: 

lconf = y/D xx T esc (6) 

From spectral modelling, we found that the mean CR density in 
the SMC is ~15% or less the CR density in the local Galaxy. 
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Assuming a similar diffusion coefficient and a potentially higher 
injection rate in the SMC, the above equations imply a mean CR 
confinement time in the SMC of ~15% the Galactic value, which 
is typically ~ 10 7 yrs, and a mean CR confinement volume size 
in the SMC at least 3 times smaller than in the MW. 

Similar conclusions can probably be drawn for the LMC. 
Discussing the possible physical reasons to this is beyond the 
scope of that paper. It may be a simple size effect like the 
surface-to- volume ratio (injection being proportional to volume 
while escape is proportional to surface), or a more sophisti- 
cated process involving the equilibrium of the various galactic 
components (CRs, gas and magnetic field resisting to gravita- 
tional attraction). It should be noted that the comparison of the 
Magellanic Clouds to the Milky Way in terms of a large-scale 
characteristic such as the CR population may be complicated by 
the fact that smaller galaxies are probably more affected by the 
stochasticity of star formation and by tidal interactions with their 
neighbours. The concept of steady-state CR population may sim- 
ply not apply in these cases. 

These rough estimates, however, hold for a mean CR den- 
sity computed under the assumption of a uniformly distributed 
CR population. As mentioned in 15.11 if CRs are poorly corre- 
lated with the gas, then the CR density can be expected to vary 
over the SMC volume and may well reach values above the lo- 
cal Galactic one. If some correlation of the CRs with supergiant 
shells exists, as mentione d in !3.3l and as was found to be the case 
in the LMC dAbdo et al.ll2009hl) . the mean CR confinement time 
and volume may well differ from the average values computed 
above. For instance, the dynamical ages of the three supergiant 
shells of the SMC are ~ 1.510 7 yrs dStanimirovic et alj|1999b . If 
CRs can be efficiently confined in the superbubbles interiors for 
this duration, the CR confinement time would be about an order 
of magnitude above the average value computed above. Such a 
scenario may be further explored observationally by looking for 
a correlation of the GeV signal with diffuse X-ray emission, the 
latter having the potential to reveal the hot gas filling the su- 
perbubble cavities. On the theoretical side, additional works are 
needed to better understand the transport of CRs in a hot, ten- 
uous and presumably highly turbulent medium and determine 
the time scale over which CRs can be retained within supergiant 
shells. 

7. Conclusion 

From the first 17 months of Fermi/LAT observations, we ob- 
tained the first detection of the SMC in high energy gamma- 
rays. The flux in the 200MeV-20GeV energy range of the data 
is (1.9 +0.2) x 1(T 8 phones- 1 , and increases to (3.7 ±1.7) 
x lO^phcm^s -1 when extrapolated to the 100MeV-500GeV 
range from our best-fit spectral model. The latter value is consis- 
tent with the upper limit derived from the CGRO/EGRET data. 

The emission is steady and from an extended source ~ 3° in 
size. It is not clearly correlated with the distribution of massive 
stars or neutral gas, nor with known pulsars or supernova rem- 
nants. However, a certain correlation of the gamma-ray emission 
with supergiant shells is observed. In the recent FermifLAT study 
of the LMC, the gamma-ray emission was observed to be poorly 
correlated with gas density but strongly correlated with tracers 
of massive star forming regions such as ionised gas, Wolf-Rayet 
stars, and supergiant shells. Better statistics from increased ex- 
posure may reveal patterns in the SMC emission and clarify the 
differences and similarities in the gamma-ray emission of the 
two Magellanic Clouds. 

If the gamma-rays from the SMC are produced by CR-ISM 



interactions, the observed flux implies an average density of CR 
nuclei in the SMC of ~15% the value measured locally in the 
MW. This value should be considered an upper limit since a sub- 
stantial fraction of the emission may be due to the population 
of high-energy pulsars of the SMC and the gas mass used for 
the calculation may be underestimated. The average density of 
CR electrons derived from radio synchrotron observations of the 
SMC is consistent with the same reduction factor, but large un- 
certainties on the radio data and on the geometry of the galaxy 
also allow higher values. From our current knowledge on the 
SMC, such a low CR density does not seem to be due to a lower 
rate of CR injection but may indicate some dependence on CR 
transport effects, such as a smaller CR confinement volume char- 
acteristic size. 
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